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Abstract 

The propagation of ultracold atomic gases through abruptly changing waveguide potentials is 
examined in the limit of non-interacting atoms. Time- independent scattering calculations of mi- 
crostructured waveguides with discontinuous changes in the transverse harmonic binding potentials 
are used to mimic waveguide perturbations and imperfections. Three basic configurations are ex- 
amined: step-like, barrier-like and well-like with waves incident in the ground mode. At low 
energies, the spectra rapidly depart from single-moded, with significant transmission and reflec- 
tion of excited modes. The high-energy limit sees 100% transmission, with the distribution of the 
transmitted modes determined simply by the overlap of the mode wave functions and interference. 
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I. INTRODUCTION 



The manipulation of ultracold matter waves can now, somewhat routinely, be performed 
above microchip or magnetized surfaces . In such experiments, the quantum nature of 
ihe dilute atomic gases dominates over the classical, enabling precision matter wave control 

One of the key requirements in using an "atom chip" to perform atom optics is the 
ability to transport atoms from one atom optical component to another. Here, we present 
calculations of wave propagation through waveguides with idealized perturbations consisting 
of sudden changes to the transverse confining potential. An increase (decrease) in the 
tranverse confining potential results in a decrease (increase) in the kinetic energy along the 
direction of wave propagation, providing effective step potentials along the waveguide. 

The present study was motivated in three ways. Firstly, recent experiments have demon- 
strated a significant fragmentation of a Bose-Einstein condensate (BEG) in a waveguide 

nn 

located close to the surfaces I5l| , attributed, at least in part, to imperfections created dur- 
ing the wire fabrication processes Q, 0, 0, [l^ ■ Secondly, to further understand some of the 
limitations to designing atom optics devices that are based on variations of the waveguide 
potentials, for example, the smoothly varying wide-narrow-wide wire geometry has been 
proposed clS 8b 

been lacking a multimode analysis of many of the simplest waveguide geometries, as it has 
been experimentally shown that introducing a perturbation in a waveguide can result in the 
transverse excitation of a BEG j^. 

To characterize the impact of transverse discontinuities, here we explore the Schrodinger 
wave mechanics of waveguides with step-like, barrier-like and well-like potentials along the 
direction of propagation. There have already been some theoretical investigations using time- 
dependent calculations of wave propagation through smooth potentials such as a bottleneck 
(step-up) and a smooth termination (extreme step-down) in the limit of non- interacting 
atoms Il^l , while non-linear (atom-atom) effects in the bottleneck- type geometry have 



also been examined 



15l |. The advantages in using abrupt potentials whilst neglect- 



ing atom-atom interactions is that simple time- independent calculations can be used to 
characterize the transmission and reflection probabilities. Under these conditions, we have 
previously investigated a circular bend , which consists of an abrupt transition from the 
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lead waveguides into the bend and at low-energies behaves like a potential well. 

Our multimode analysis, restricted as it is to the the linear regime, provides a baseline for 
comparison of BEG propagation through quasi-one- dimensional (1-D) waveguides including 
the transverse degrees of freedom. For example, previous time-independent studies have 
investigated non-linear wave propagation through shallow-well, step and gaussian shaped 1-D 



potentials |l2 , llSl ll9|, |20| . Such simple waveguide potentials could be generated by modifying 



the transverse confinement, where knowledge of the transverse excitation probabilities, in 
the abrupt and linear limits, should be useful. 

For ground mode matter waves propagating at low energies through the various pertur- 
bations, the present results show that the spectra rapidly depart from single-moded, with 
significant transmission and reflection of excited modes. The high-energy limit sees 100% 
transmission, and we present a simple model to determine the distribution of the trans- 
mitted modes that combines the overlap of the mode wave functions with the multi-path 
interference of the modes. 



II. DETAILS OF THE CALCULATIONS 

There are a number of atom chip wire configurations that can create waveguides 
A bu. we fonow t.eo.e^cal a.a. adopted ,„ Re.. U Q- Tha. we a.u.e 



that the waveguides consist of an idealized trapping potential that is quadratic near the 
minimum and operate at low matter-wave densities such that atom-atom interactions can 
be neglected. Furthermore, so that the waveguide potentials reduce to an effective 2D 
problem, the waveguides are assumed to be created by multiple wire configurations with 
abrupt changes in the spacing between the wires, such that the height and transverse position 
of the potential minima remains constant along the waveguide. The out-of-plane quantum 
number is then a conserved quantity. 

We consider simple harmonic oscillator (SHO) confining potentials, which, for barrier-like 
or well-like effective potentials, are given by 



V{x,z) 



1 2 2 



Zo < Z < Zi , (1) 



k 2 



^mu'^x^, z > zi . 
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The barrier-like potential has oob > oOa, the well-like potential has oOa > uji,. The step- like 
potential consists of only one change in frequency. 

Oscillator units are used throughout this paper, where energies are in units of hu, while 
lengths are given in units of /9 = ^JhJmhJ. An example barrier-like potential is shown in 
Fig. m^a), where the reference frequency is tUa = 1 and ojh = 2. 

(b) 



FIG. 1: (a) Potential energy surface of a barrier-like waveguide with cJq = 1 and w;, = 2 and barrier 
length / = zi — zq. (b) Energy levels of the leads and barrier transverse SHO potential along z. 
The solid lines at zq and zi should be vertical, but instead are drawn on an angle to highlight 
the lead-barrier-lead transition. The dot-dashed lines correspond to the parity-forbidden levels 
(assuming an even incoming mode). All the energies and coordinates are given in terms of oscillator 
units for the leads. 

The corresponding energy levels of Eq. (P) are shown in Fig. ^b). These energy levels 
behave as effective potentials for the longitudinal motion since we expand the total wavefunc- 
tion in each region on the transverse oscillator states. In this model, all coupling between 
modes occurs through the matching between regions. The potentials in Eq. are sym- 
metric in X so that parity in x is conserved, simplifying the present analysis and discussion 
considerably [the dot-dashed lines in Fig. C^b) are not coupled to the solid lines]. In ex- 
periments, imperfections would as likely be off-center, resulting in populating of all modes 
(as was possible in our previous study of the circular bend The fundamental physics, 

however, remains much the same, so we chose to adopt parity-conserving perturbations with 
the incoming waves restricted to the ground (even) oscillator mode. 

To perform the time-independent scattering calculations we initially adopted the transfer 
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matrix method 23l | , al though most of the calculations reported in this paper use the interface 
matching method [l6|, |23|. The two methods are similar, however, and since the transfer 
matrix approach facilitates the discussion of our results, we outline it here. The extension 
of the transfer matrix method from one-dimension to include transverse degrees of freedom 



is trivial 



2J, |25[, so only a short summary is given here as it applies to a barrier/well-like 
geometry of length {zi — zq) = I. 

Firstly, the time- independent wavefunction is expanded on transverse SHO states, ipn{x) 
for frequencies uja, and Xm{x) for Ub. 

n 

^n(x, z)=Y, Xm{x) [c^e^-^'-'°^ + d^e-'^-^'-^'^] z, < z < z^ , (2) 

m 

v|/in(x, z)=Y, '^nix) [gne'""^'-''^ + Ke-'^-^'-'^^] z, < z . 

n 

where the momenta are = \/2E — ujai'^n + 1) and Km = ^/'2E — Uhi^m + 1). Matching 
the wavefunctions and their first derivatives across each interface and then projecting out 
the modes gives the following sets of equations: 



: = " : J ! ■ (3) 



The matrix elements of each submatrix are 

411 -422 _ 1 A I ^ \/n 





412 _ /i2i _ 1 A _ ^ \/n 

^nm ~ ^nm ~ 2 V k ' 



nm } 



(4) 



^mn ~ ^ ^mn ~ 2\ k ) 

0l2 _ -2iKmln21 — \ ( T _ p-it^ml ri 

^mn ~ ^ ^mn ~ 2 V k ) '^nm , 

using the notation Onm = J '^n{x)Xm{x) dx. While generating functions are known for the 
overlaps of SHO-functions with different frequencies Q^El (such overlaps are also found 

28|), we performed the 



in calculations of transitions between molecular vibrational modes 
transverse integrations numerically using a B-spline basis. 

The scattering solution is obtained by constructing the transfer matrix Q = AB, which 
relates one lead's coefficients, a and b, to the other's, g and h. Given that the wave is 



restricted to incoming from z < Zq in the ground mode {qq — 1, a„.>o = 0, and h — 0), 
the hnear equation a — Q^^g is solved, and then b — Q^^g. The transmission and reflection 
probabihties for each mode are then given by: 

\n I'^h \h I'^h 

Tn/ = I and R„ = ' . 5) 

There are convergence difficulties with the transfer matrix approach, the demonstration 
and discussion of which is mostly relegated to the appendix. In brief, the problems are related 
to the slow decay of the SHO overlaps Oom with m, which requires the inclusion of strongly 
closed channels, leading to the appearance of large exponentials in the transfer matrix. While 
the transfer matrix method eliminates the need to find intermediate coefficients (ic. c and d), 
it was generally found to be unstable beyond the smallest of the perturbations considered in 
this paper. The interface matching method, in contrast, explicitly solves for the intermediate 
coefficients and is able to include enough closed channels to ensure near machine-precision 
convergence for the range of geometries and energies given in this paper. The transfer matrix 
method works well, however, for step-like geometries with a single interface (ie. Q = A) as 
there are no exponentials in the matrix elements, and furthermore, all the coefficients are 
explicitly solved for. 



III. RESULTS 



The multimoded transmission and reflection probabilities for three basic geometries are 
given in this section: the step-like potential, the barrier-like potential and the well-like 
potential. The calculations for the step-like potentials use the transfer matrix method, 
while interface matching is used for the barrier-like and the well-like potentials. Before 
presenting these calculations, it is instructive to discuss how the interface overlaps Oom 
scale with frequency. 



A. Interface overlaps 

The Oom dependence on the frequency ratio ooh/ooa is discussed here since it strongly 
influences the amount of mode excitation caused by the different geometries. Using the 
orthonormality of the SHO eigenstates, along with the recursion relations of the Hermite 
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polynomials, it can be shown that the overlap integrals of the uja = 1 ground mode with the 
Ub even modes reduce to the particularly simple form: 



Om 



(6) 



{uJb + 1) m\\ \uJb + 1 

for all m = 0, 2, 4, 6 . . .. Due to symmetry, Oom = for m = 1, 3, 5 . . .. Eq. (jH} was also 
obtained by Aslangul j^l- 

The dependence of the first three overlaps (m = 0, 2, 4) on frequency is shown in Fig.|2fa). 
At UJb = 1, there is a perfect waveguide match, and we must have Oqo = 1 and Oom = for 
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FIG. 2: SHO wavefunction overlaps |OomP between the ground mode with fixed frequency uJa = 1, 
and another mode with variable frequency tob- (a) Shows the m = 0, 2, 4 modes as a function of 
iOb- The arrows labelled by uJm indicate the maxima, (b) Shows the overlaps as a function of m for 
four frequencies: tOb = 1.5 (squares), tob = 2 (filled-triangle) , tob = 5 (hollow-triangle) and tob = 20 
(circles). The lines between dots in (b) are added to guide the eye. The frequencies shown in both 
(a) and (b) are given in oscillator units relative to Wa- 
rn > 0. As UJb increases, Oqo monotonically decreases towards zero. At the same time, the 
overlap with each excited mode increases and reaches a maximum when the characteristic 
equation, — 2(2m + l)^;^ + 1 = 0, is satisfied. As a function of frequency, these maxima 
occur at lu^ = (2m + 1) + a/ (2m + 1)^ — 1, which, for m = 2, is 0^2 = 9.8989 and for m = 4, 
is 004 = 17.944. Comparing the width (a/ {x'^)) of the SHO functions at these frequencies 
against the ground mode's reveals that ^y {x^iuj^)) ~ a/1/4 = a/ (xq)/2. This connection, 
while natural, is not particularly illuminating and is not pursued here any further. Past the 
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maxima, the overlaps slowly decrease as to zero. Due to symmetry of the ratio Uh/uJa, 

a second maxima of |OomP also exists at l/ujm- 

The slow decay of the overlaps with the higher modes can be seen in Fig. Efb). This 
implies that, even at low energies, many closed channels must be included in the following 
calculations to achieve computational convergence. 



B. Step-like waveguide potential 

The transmission and reflection probabilities, T„^ and R„^. , of ground mode plane- waves 
traversing four step-like waveguides are shown in Fig. IHl At incident energies below the 
lowest reflection threshold {E < 2.5, the lowest excited mode energy), the system behaves 
like the familiar 1-D step potential. Ground state transmission, Tq, remains the dominant 
channel across the range of energies shown in Fig. IHl^a) (for the range of Ub examined here), 
although excited mode transmission can also be seen in Fig. Efa) as each mode opens. 
Significant reflection is seen into the ground nj = mode in Fig. IHfb), which rapidly drops 
off from threshold. As each reflection threshold opens, the reflection into the excited modes 
Uf > 0, seen in Fig. Ofc), firstly increases then is seen to experience an overall decrease. All 
of which are consistent with the Wigner threshold laws for multichannel systems js^. 

For incidence in the ground channel, the total transmission approaches 100% in the high 
energy limit. In this limit the transmission probability is given simply by 

T„,(E^oo) = |0o„/ , (7) 

with the SHO overlaps of Eq. Such projections were introduced as part of the waveg- 
uide calculations of Jaaskelainen and Stenholm Q,!^^, in which the transmission excitation 
probabilities generated by both bottleneck and split-potential waveguides were briefly dis- 
cussed as the potentials tended towards abrupt. Similar multimode projections have also 

)een theoretically examined during expansion of a BEG from a microtrap into a waveguide 

15^ , with an emphasis on the effects of atom-atom interactions. 

Equation (jZj) can be seen as the high-energy limit of the matrix elements of A, given by 
Eq. (Q. At energies high compared to the step height, the momenta of the lowest few modes 
are approximately the same on both sides of the step, fi;„^ ^ k^. Given that the overlaps limit 
the number of channels involved, the matrix elements that dominate the Q = A transfer 
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FIG. 3: Multimoded scattering probabilities of step-like potentials which consist of a single abrupt 
change in the transverse confinement potential from = 1 to a;;, = 1.5,2,3, and 5. The incoming 
waves are in the ground state = 0. (a) Shows Tq and ^ T„^.>o, the transmission probabilities into 
the ground mode and the sum of the transmission probabilities into the excited modes, respectively, 
(b) Gives Rq, the reflection probabilities into the ground mode, while (c) gives X]Rnj>0; the sum 
of the reflection probabilities into the excited modes. The total energy E is given in oscillator units 
relative to uJa- The crosses in (a) and (b) are the analytic transmission and reflection probabilities 
for a 1-D step potential for the case of uj^ = 5 (ie. of height Vq = |2.5 — 0.5| with a background 
potential V = 0.5 added to correct the reflection ground mode threshold). 

matrix are then ^on^ — ^o\f ~ Oouf, while ^o^f ~ ^On^ ~ 0- For an incoming wave in the 
ground mode with = 1, a„.>o = 0, and d = 0, then the outgoing waves have ~ Oon^ 
while 6 = 0. 

To more clearly show this limit, the results of Fig. El are replotted in Fig. H] on an energy 
axis scaled by Uf, instead of Ua so that the transmission channels for each waveguide open at 
the same scaled energy. At high energies, the transmission probabilities shown in Fig. Ufa) 
increase towards asymptotes of To{E — > oo) = |OooP, in agreement with Eq. 0. For the 
four different waveguides shown here, the limits are To{E —>■ oo) = 0.979796, 0.942809, 
0.866026, and 0.745356. Figures lU^b) and (c) show that at the onset of each transmitted 
mode (ie. at E/Ei, = 5,9 etc.) the transmission probability into that mode increases, 
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FIG. 4: The scattering probabilities of the same four step-Hke potentials as Fig. (w;, = 1.5,2,3, 
and 5) plotted here as a function of E/Eh, where Ef^ = Wf,/2. For each waveguide, (a) shows 
the transmission probabilities of the ground mode, (b) the sum of the excited mode transmission 
probabilities X]'^n/>0) (c) the total reflection probabilities ^R„^. The arrows at E/Ef^ = 15 
correspond to the |OooP interface overlaps. 

taking flux from reflection. The mismatch in mode wavefunctions for uJb = 5, for example, is 
particularly severe, with T„^=o,2,4,6(^ ^ oo) = 0.745356, 0.165635, 0.055212, and 0.020449 
[see Fig. El^b)]. In this case, these four transmission modes must be open before the high 
energy limit T„^, = 1, as per Eq. ((7j)] is reached to within 2%. 

Whether the wave is incident from the left or the right, these results apply. The ground 
mode transmission To(-E') is absolutely identical as a function of total energy E for both the 
waveguide constriction (ie. step-up from = 1 to cuf, > 1), and the waveguide expansion (ie. 
step-down from Ua > 1 to Ub = 1) This results also holds for the familiar 1-D step potential, 
and, although we do not show it here, the transmission and reflection mode mixing conspires 
to ensure this is also the case in the multichannel system. The mode mixing as a function 
of E is different for either geometry, however, since for the step-up case there can be many 
reflection channels open at the lowest transmission threshold, E = Eb, while for the step- 
down, there can be many transmission channels open dX E = Eb- At high energies, neither 
of the step geometries generates reflection, and for incidence in the ground mode the limits 
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from Eq. (|7j) apply. 



C. Barrier-like waveguide potential 

To demonstrate the characteristics of a single barrier-like potential, we consider the case 
shown in Fig. [H for a fixed length / and frequencies that change from Ua = I to Ub > uja 
and back to Ua = 1- We also present the high-energy characteristics for scattering from this 
potential. 

The transmission and reflection probabilities for four waveguide constrictions ujfy = 
1.5,2,3,5 and length I = 10 are shown in Fig. El Resonances appear in all of the spec- 
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FIG. 5: Multimoded scattering probabilities of four barrier-like potentials with Uh = 1.5,2,3, and 
5 for a fixed length / = 10 as a function of E (in oscillator units relative to uJa)- The legend for 
the different frequencies is the same as Figs. (jSJ and Q. The organisation of the probabilities is 
also the same as Fig. The crosses in (a) and (b) are the analytic transmission and reflection 
probabilities for a 1-D barrier potential for the case of u;^ = 3 (ie. of height Vq = |1.5 — 0.5| with a 
background potential V = 0.5). 

tra at low energies due to the wavelength matching condition nX/2 ^ /. The 1-D analytic 
results for cj^ = 3 in Fig. Efa) highlight the transmission resonances as the result of low- 
energy ground- mode propagation over the barrier. As soon as the Ua = 1, n = 2 mode 
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opens at E = 2.5, however, multichannel physics takes over (n = 1 is not allowed due to 
symmetry). Above E = 2.5, there is significant excited mode transmission [Fig. |Sfa)], while 
the amount of refiection into the ground mode [Fig. E^b)] and the excited modes [Fig. Etc)] 
is significant across the energy range. 

To observe the high-energy limit, the results for the cuj, = 5 barrier of Fig. El were extended 
to higher energies, and are shown in Fig. El with the energy scaled by Ei, = 2.5. Figure Efb) 



A 
C 
I- 

o 

c 

I- 



o 

A 

CC 
o 

o 
DC 

FIG. 6: The scattering probabilities of the iVh = 5, I = 10, barrier-like potential seen in Fig. El 
plotted here as a function of E/ E^, where E^ = 2.5. (a) Shows the individual transmission proba- 
bilities of the Uf = 0, 2, 4 modes alongside the sum of the transmission probabilities for the Uf > 4 
modes, (b) Shows the ground-mode reflection probabilities as well as the sum of the excited mode 
reflection probabilities. The crosses in (a) correspond to the To(-E' — > oo) two-mode interference 
model of Eq. The first peak of Tq reaches nearly up to 1, which is not shown due to the limited 
energy resolution of the figure. 



shows that while there is more structure in the refiection probabilities than for the step-like 
potential in Fig. 01^c), the total amount of refiection still tends towards zero as more barrier 
modes become open at E/Eb = 5, 9, 13, . . .. 

To obtain an expression analogous to Eq. ((7|), we must take into account the fact that 
the transmitted waves going through a barrier experience at least two interface projections 
as in Eq. The phases accumulated while propagating the length / of the barrier must 
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also be included in such a prescription, which suggests that at high energies (relative to the 
barrier height) 

To(^ oo) = I ^ e^(''™-'^'')'OomO™of , (8) 

m 

where m is only summed over the propagating barrier modes. The crosses shown in Fig. IHta) 

demonstrate that the two-mode version of this model does a remarkable job in describing 

the transmission probability above the m = 2 barrier threshold (at E/Eb = 5). Depending 

on the phase differences, the modes that are excited at the first interface can be converted 

back to the ground mode by the second interface. We also noted this behavior in circular 

waveguide bends where the amount of excitation could be suppressed by changing the 

angle swept out by the bend to the point where the accumulated phase difference between 

the n = and n = 1 modes was a multiple of vr. A similar design consideration could 

I I 

perhaps be useful for atom optics devices such as the quantum point contacts llj, where 
any unavoidable - yet unwanted - mode excitations could be minimised by varying the 
length between the changes in waveguide potentials. 



D. Well-like waveguide potential 

The transmission and reflection probabilities for well-like waveguides due to a potential 
bulge > ujb) are the focus in this last section. The scattering behavior for a well is 
complicated by the presence of bound states which translate into the presence of Feshbach 
resonances in a multichannel problem. Much of the resonance physics seen here has been 



extensively discussed as part of our studies of the circular waveguide bend 
propagation thresholds in a bend lie slightly lower than the connecting leads 



32 



m 



3l|. The 

result- 



ing in very weakly bound states and energetically narrow resonances. The present well-like 
waveguides can provide extreme differences between the lead and bulge energy thresholds, 
and thus the possibility of multiple narrow resonances located below the thresholds. 

Two bulges are considered here, the first from uja = 1-5 to Ub = 1 and back to uJa = 1.5, 
and the second with Ua = 2 to Ub = I and back to Ua = 2. In both cases, we choose ujb = I 
to be the reference oscillator frequency, to simplify the comparison with the barriers in the 
previous section. The transmission and reflection probabilities of both of these potentials 
with well length / = 10 are shown in Fig. [7| 

Multiple 100% reflection resonances exist at energies below the first excited lead mode 
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FIG. 7: Multimoded scattering probabilities of two well- like potentials with coa = 1-5 and 2 for a 
fixed length / = 10 as a function of E (in oscillator units relative to ujh). The organisation of the 
probabilities is the same as Fig. 01 The crosses in (a) correspond to the To(-E — > oo) two- mode 
interference model of Eq. @ for uJa = 2. 

energy, ie. below E = 3.75 for uJa = 1-5 and E = 5 ioT uJa = 2. As was noted in our 
circular bend studies the reflection resonances below the second excited mode (ie. 
below E = 6.75 for Ua = 1.5, and at -E = 9 for Ua = 2) do not result in complete reflection 
due to the reduced coupling between the ground and second excited mode and due to the 
existence of alternate pathways to reflection. 

The high-energy transmission probability asymptotes for the well-like potential are again 
given by Eq. (jH)), although the energy of the first excited mode within the well means that the 
two-mode model starts at E = 2.5 (in other words, at a lower energy than for the equivalent 
barrier-like potential). The two- mode model is shown in Fig. [3^a) for the Ua = 2 well, and is 
seen to be a bad approximation near the E = 2.5 threshold due to the significant reflection 
there. The two-mode model generally provides a reasonable approximation at higher energies 
once the total reflection probability has dropped below the total excited mode transmission 
probability, and also when the many narrow reflection resonances no longer play a role. 

Waveguides with more extreme discontinuities such as = 3 and 5 were also explored. 
They exhibit so many resonances across the range of energies shown in Fig. however, that 
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the transmission and reflection probabilities essentially become a dense series of closely- 
spaced vertical spikes. At energies above the second excited mode threshold {E = 4.5), 
there are not so many resonances, although there remains significant reflection probability. 
As an example, the Ua = 5, I = 10, well-like waveguide has a total reflection probability at 
E = 37.5 {E/Ea = 15) that has only dropped down to about 0.2. This reflection probability 
can be compared with the Ub = 5, I = 10 barrier-like geometry in Fig. IHl^b) at E = 37.5 
[E/Eb = 15), where the total reflection probability was only about 0.02. In other words, 
reflections play a far more dominant role for the well-like potentials than the barrier-like. 



IV. SUMMARY 



Using time-independent solutions of the Schrodinger equation, we have explored the 
propagation of dilute, ultracold atomic gases through abruptly changing waveguide poten- 
tials. Previous studies have discussed the conditions for "adiabatic" waveguide propagation 
through microstructures (eg. Refs. In contrast, the interest here was on 

the consequences of sudden potential variations for mode excitation with a view towards 
modelling waveguide imperfections, examining the effects of using abrupt potentials in atom 
optical devices, and simply to explore the behavior of some simple geometries. 

Three idealized geometries with changes in the transverse guiding potential have been 
the focus of study: step-like, barrier-like and well-like. The low-energy behavior of all the 
geometries departed from single-moded, with the exception of the mildest perturbations at 
energies below the lowest excitation threshold. Significant generation of both transmission 
and reflection excited modes was caused by the mismatch of the modes at the interfaces 
between the waveguide sections. The strong coupling to excited modes is due to the signifi- 
cant overlap of the ground SHO function of one frequency with the excited SHO functions 
with a different frequency. Certainly, care should be taken during wire fabrication of atom 
optical elements to ensure that any deliberate (or not) changes in a waveguides transverse 
frequency are not abrupt. 

High energy wave propagation through abrupt potentials amounts to 100% wave trans- 
mission via projections across each interface, along with multimode interference. The present 
SHO-based waveguides behave somewhat differently than the hard-walled models for ballis- 
tic electron propagation through waveguides with a wide-narrow junction (see Ref. Q 

, and 
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the references thereafter that cite it). In that case an impedence mismatch occurs, where 
there is always some amount of reflection at high energies due to the fact that the narrow 
guide modes can never represent the wide guide modes over the whole width of the lead. 
Although, to the best of our knowledge, calculations for such electron waveguides have not 
discussed high-energy /abrupt potential transmission and reflection limits. 

A further condition was suggested for high energy transmission through microstructures 
with multiple interfaces to account for the interference between modes. It was shown that 
a simple two-mode model can give a reasonable approximation to the amount of ground 
mode transmission, and provides an additional consideration for the design of "atom-chip" 
waveguides to control single-moded wave-transmission through potentials that generate mul- 
timoded excitations. 
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APPENDIX A: CONVERGENCE OF TRANSFER MATRIX CALCULATIONS 



Our transfer matrix program was validated by reproducin g th e single-mode calculations of 
electron propagation through a linear array of 1-D potentials j36[ , and secondly by comparing 
with multi-moded 2-D results from an interface matching program IgI ]. 

The multimode transfer matrix method, however, has numerical convergence difficulties 
as the number of modes included in a calculation increases. This behavior can be seen in 
Table HI for scattering with different incident energies off a barrier-like potential with 
ujii/ uOa = 2 oi length / = 5. The unitarity of the transfer matrix results deviate significantly 
from 'Ylif'^uf + = 1; even as is increased. The interface matching calculations are 
shown for comparison, and barely suffer from the same problems. 

In an attempt to rescue the transfer matrix approach, a second method for constructing 
the transfer matrix was investigated. This second attempt involved inverting a matrix C, 

-1 / \ 

(Al) 
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TABLE I: The convergence of transfer matrix and interface matching calculations for a barrier-like 
potential of length / = 5 with uo^l^a = 2 The number of channels included in each calculation is 
given by N. The ground mode transmission probabilities and unitarity (given by 1— ^j(T„^+R„^)) 
are given for three energies: E = 0.75, where only barrier tunnelling is energetically allowed, and 



E = 6, where three lead modes and two barrier modes are open. 







E = 0.75 




E = Q 




AT 


lo 


Unitarity 


lo 




Unitarity 






Analytic (1-D barrier) 








1 


0.00339154 




0.999976 










Transfer Matrix Q = AB 








1 


0.00381548 


-4.2394-4 


1.124973 




-1.2500-1 


2 


0.00267029 


3.0328-5 


0.772832 




-5.2616-2 


4 


0.00274164 


2.8303-'^ 


0.804769 




8.1844-4 


8 


0.00274223 


4.1060-'^ 


0.807601 




-3.4533-* 


16 


0.00264476 


-9.4352-3 


0.849707 




-3.4386-2 






Trantifpr M^itrix O — AC!^^ 

_L 1 CLiii!!)is_.i IVXOj Ui l-V Kc^ — 








1 


0.00339154 


3.3307-1'^ 


0.999976 




2.0912-16 


2 


0.00270062 


-1.4129-12 


0.795176 




3.3826-2 


4 


0.00274192 


1.7813-1° 


0.805169 




-3.1198-3 


8 


0.00274223 


-4.5390-® 


0.807546 




-1.5052-5 


16 


0.00277009 


4.0555-3 


0.916579 




-5.6275-2 






Interface Matching 








1 


0.00339154 


2.2204-16 


0.999976 







2 


0.00270062 


1.1102-i'3 


0.795176 




3.3826-2 


4 


0.00274192 


3.3307-16 


0.805169 




-3.1198-3 


8 


0.00274223 





0.807605 




-3.8357-6 


16 


0.00274223 


4.4409-16 


0.807607 




-7.6739-13 


32 


0.00274223 


1.1102-16 


0.807607 




2.2205-16 



with matrix elements, C^^ = = 1(1 + Km/kn)Onm, and C^^ = = 1(1- K,m/kn)Onm 
while the matrix elements in A of Eq. Q were modified to include the e^*''™' terms. The 
numerical results of solving Q = AC^^ are also shown in Table HI and while unitarity results 
for the single-channel calculation, it is as unstable, if not worse that Q = AB. 

For the geometries considered in this paper the transfer matrix approach was not able to 
adequately converge over the required range of energies. The exception to this was for the 
step-like geometry, where for the most extreme calculations in this paper (Fig. |31 uJb/oJa = 5 
and E = 37.5) = 128 modes could easily be included (since the transfer matrix contains 
no exponential terms, and all coefficients are explicitly solved for). 
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We also find that, since the Onm overlaps are energy independent, obtaining convergence 
is much easier for the present interface matching calculations than for previous circular 



16,n7||. For the most extreme barrier- 



waveguide bend calculations using the same method 
like calculation (Fig. IHl ujb/^a = 5, / = 10 and E = 37.5), the interface matching method 
achieved unitarity to better than 10^^° with = 96 modes. 



Some of the transfer matrix convergence issues were passingly mentioned by Wu et. al j2J| , 
where the transfer matrix method was initially applied to a system of multiple circular bends 
(the transfer matrix is simply built as a number of matrix products). In the end, however, the 
calculations shown in Ref. 2^1 were performed using multiple interface matching. We labor 
on this point since no such unitarity problems were noted in the multichannel transfer matrix 
calculations of Ref. |2^], where the examples given were for 2 or 3 channel calculations with 
delta potentials. On the other hand, the advantage of the transfer matrix approach is that, 
in the single-mode approximation, one can easily obtain the transmission characteristics of 
both finite and infinite length periodic potentials 
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